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Abetract 

Thia  work  praeenta  an  alternative  to  current  integral 
approaches  that  do  not  accurately  predict  noise  from  a  de¬ 
localised  rotor.  The  two-dimensional  model  of  this  study  is 
essentially  an  airfoil  shape  mounted  on  the  side  of  cylin¬ 
der  which  spins  at  transonic  speeds  between  two  end- 
walls.  A  Unite-difference  conservative  formulation  of  the 
two-dimensional,  transonic  full- potential  agnations  solves 
the  flow  field  sway  from  the  body  directly  and  hence  pre¬ 
dicts  the  beginning  of  delocalisation  and  the  associated 
wave  propagation.  The  potential  method  has  an  inherent 
stability  limit  at  a  free  stream  velocity  of  Mach  1.4  and 
shows  a  troubling  sensitivity  to  the  choice  of  outer  bound¬ 
ary  conditions.  Methods  for  correctly  handling  these  two 
problems  are  shown.  A  two-dimensional  experiment  that 
simulates  the  conditions  assumed  by  the  computer  coda  is 
described.  The  computer  prediction  of  the  flow  expected 
in  the  experiment  is  presented  as  are  some  limited  exper¬ 
imental  results.  Results  taken  from  a  three-dimensional 
hovering  rotor  code  show  how  the  current  two^Umenshmal 
results  correlate  to  those  from  an  actual  rotor.  Results  for 
low-aspect-ratio  rotors  show  a  higher  delocalisation  Mach 
number  for  decreasing  aspect  ratio  and  a  greater  depen¬ 
dence  on  thickness.  ' — 

Introduction 

High-speed  impulsive  (HSI)  noise  is  emanated  whan 
a  transonic  rotor  creates  sack  a  strong  shock  that  the 
shock  leaves  the  blade  and  propagates  into  the  free  stream 
This  free-stream  shock,  or  delocalised  shock,  creates  sev¬ 
eral  problems  for  integral  methods.  Earlier  research  has 
been  applied  to  this  problem  but  with  little  success.1-* 
More  recent  works  have  fared  better4-*  with  this  ‘delo¬ 
calization*  problem,  as  Caradonna  and  Isomc  termed  this 
spreading  of  the  shock.  Photographs  of  such  delocalised 
shocks  radiating  from  model  propellers  were  produced  by 
Hilton10  in  1938  and  more  recently  by  Tangier11  in  1977 
for  a  model  helicopter  rotor  in  a  forward  flight  condition. 
A  similar  phenomenon  is  seen  in  hovering  helicopter  ro¬ 
tor  tests  at  high  tip-Msch-numbers.  Hot-wire  studies  from 
Schmitz  and  Yu1*  experimentally  demonstrate  the  extant 
of  the  supersonic  zones  relative  to  the  rotor  on  and  around 
the  tip.  Figure  1  shows  results  from  their  study  along  with 
the  associated  changes  in  sound  signatures  as  the  tip-shock 
delocalises.  These  results  are  from  a  hovering  rotor  model 
with  a  NACA  0012  airfoil  and  aspect  ratio  of  13.7. 

One  of  the  moat  familiar  integral  formulations  for  solv¬ 
ing  acoustics  problems  was  presated  by  Ffowcs- Williams 
and  Hawkings:1* 
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Applications  using  the  linear  terms  of  this  equation  work 
well  except  wha  a  large  shock  appears  in  the  flow.14  Eva 
using  the  full  nonlinear  portion  of  the  Ffowcs- Williams  ad 
Hawkings  (FWH)  equation  gives  incorrect  trends  wha  a 
rotor  has  a  delocalised  tip  shock.1*  Schmitz  ad  Yu14  use 
the  first,  or  monopole,  term  of  the  FWH  equation  in  their 
linear  analysis  of  a  hovering  rotor.  This  dipole  term  has 
almost  no  effect  on  the  predicted  pressure  pulse  wha  com¬ 
pared  with  the  monopole  or  thickness  term  prediction.14*1* 
The  third  ud  last  term  of  the  FWH  equation  is  the 
quadrupoie  term,  which  has  no  simple  physical  analogue 
like  the  other  two  terms.  Quadrupoie  terms  describe  the 
fluid  distortions  ud  pressures  over  a  volume  near  the  blade 
tip  ad  are  quite  necessary  to  accurate  tracking  of  a  dis¬ 
turbance  as  it  moves  off  the  blade.  This  last  term  is  ex¬ 
tremely  difficult  to  integrate  or  model,  but  its  inclusion  has 
baa  shown  to  improve  the  pulse  predictions  significantly. 
Yu  et  al.18  completed  a  complex  nonlinear  quadrupoie  so¬ 
lution  that  greatly  improvae  the  solutions  but  still  ha 
difficulty  with  delocalised  cases.17  Perhaps  a  differential 
method  would  work  better  in  such  cases. 

There  are  several  three-dimensional  codes  that  model 
the  transonic  rotor-tip  flow  field  using  flnite-diflerence  (FD) 
formulations, ia,w  but  these  codas  only  predict  the  flow  near 
the  rotor  tip  ud  tend  to  damp  any  waves  that  spread  off 
the  blade  surface.  Kirchhoff  methods  map  a  linear  pressure 
wave  to  some  other  physical  point*0  wha  provided  with 
a  accurate  starting  wave.  However,  there  ie  no  direct  way 
to  connect  the  surface-pressure  solution  with  the  praaeure 
field  off  the  blade. 

There  has  baa  one  research  effort  that  couples  a 
two-dimensional  FD  solution  near  some  arbitrary  dtetur- 
bance  to  u  acoustic-wave  solver  away  from  the  blade. 
Rutherford*1  praeatad  a  small-disturbance,  noocooservs- 
thre,  two-dimensional  shock-fitting  potential  code  for  the 
innermost  portion  and  a  MacCormack  solution  of  the  dis¬ 
cretized  wave  equation  which  uses  initial  boundary  condi¬ 
tions  defined  by  the  potential  code. 

The  work  presented  hare  develops  experimental  ud 
numerical  techniques  for  the  direct  prediction  of  noise  prop¬ 
agation  from  a  two-dimensional  rotating  cylinder  with  a 
potential  FD  method.  This  two-dimensional  approach  al¬ 
lows  investigation  of  such  problems  as  stability  limits  and 
the  selection  of  proper  boundary  conditions,  problems  that 
arise  in  potential  codes  used  in  supersonic  sones. 

Previous  studies  by  Schmits  and  Yu,1*  which  show 
that  only  the  last  10X  of  a  blade  determines  most  of  the 
HSI  noiee  levels  ud  that  lift  has  little  effect,  indicate  the 
acceptability  of  the  two-dimensional  approach.  Figure  2 
shows  how  the  two-dimensional  cylinder  model  correlates 
to  a  generic  rotor  blade.  Such  a  two-dimensional  exper¬ 
iment  is  described  ud  conducted  to  verify  the  simplified 
model  ud  computer  code.  Figure  3  is  a  sketch  of  this  phys¬ 
ical  test.  Boundary-layer  problems  found  while  running  the 
tests  mandated  a  correction  to  the  inviscid  computer  code. 
Limited  results  from  this  experiment  ud  the  computer  so¬ 
lutions  are  presated  to  validate  the  code  ud  to  form  a 
basis  for  extrapolations  to  real  helicopter  problems.  Some 
surprising  trade  appear  in  the  predicted  results  for  the 
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short  rotors  typical  of  a  propeller  or  tilt- rotor- type  appli¬ 
cation. 


The  usual  unsteady  continuity  and  density  equations 
in  cylindrical  coordinates  provide  the  bask  for  the  finite- 
difference  code  of  this  work.  After  transformation  to 
a  routing  coordinate  system,  a  computational  variables 
transformation33  yields  the  following  nondimensional  equa¬ 
tions  for  programming: 

(pU)t  +  (pV)r-  0  (2) 

+  (3) 

where  the  velocities  U,  Ugo,  and  V  are  defined  as 
(Y  +  Ads),  Y,  and  [Y4r)r,  respectively,  and  where  p  k  the 
fluid  density,  q  k  the  ratio  of  specific  heats,  a  k  the  free- 
stream  speed  of  sound,  d  k  the  usual  flow  potential,  r  k  the 
radial  dimension  and  distance,  and  $  k  the  circumferential 
dimension. 

The  density  p  and  the  velocity  components  £d*  and 
dr  are  nondimenshmalised  by  ambient  density  ana  by  the 
cylinder  surface  velocity,  OR,  respectively,  where  H  k  the 
angular  velocity  and  R  k  the  ramus  of  the  cylinder.  The 
radial  dktance  r  k  nondimenskmalised  by  chord  length  c; 
the  circumferential  variable  t  k  nondimensionalked  by  the 
sector  width  occupied  by  the  airfoil  on  the  cylinder  edge. 
Thk  sector  width  k  also  given  by  the  inverse  of  the  aspect 
ratio,  At,  which  k  similar  to  the  aspect  ratio  of  a  rotor  blade 
and  k  defined  as  R/e  in  previously  defined  variables.  An 
additional  factor,  Y,  k  introduced  to  convert  the  scale  of 
the  radial  dktance  variable  to  cylinder  radii  lengths.  Thk 
factor  Y  k  defined  as  r/R,  which  appears  as  the  radial 
dktance  in  the  equations  below.  Finally,  »  (IR/a  k 
the  cylinder  Mach  number. 

Further  transformations  for  computational  consid¬ 
erations  follow  conventional  procedures  for  conservation 
forms:"” 

(*).♦(*).-  ^  « 
p  *  [l "  ^ +  2B4'+' +  c*?)]  H  1  {5) 

The  quantities  U  and  V,  called  the  contravariant  veloc¬ 
ity  components  along  the  (  and  q  directions,  respectively, 
are  defined  below.  Metrics  A,  B,  C,  D ,  and  E  and  the 
Jacobian,  J,  are  groupe  of  various  mapping  quantities  that 
arise  during  transformation;  they  are  defined  as 

A  »  +  y(r 

B  -  y(tn*  +  V(r1r 

c  -  ynl  +  Yr,i  («) 


j  -  (sOr  -  frU# 


A  sue  case  ire- line  over-relaxation  (SLOR)  method 
solves  the  differenced  equations  outlined  below.  An  arti¬ 
ficial  density  scheme  builds  in  proper  tones  of  dependence 
and  allows  the  use  of  simple  central-difference  operators  in 
all  regions.  Second-order-  accurate  formulas  determine  the 


flux  terms  fU  and  pV  at  midpoint  locations.  The  con¬ 
tinuity  equation  k  differenced  to  males  the  total  operator 
centrally  differenced  about  a  node  point: 


*(?)♦*(?)- 


where  6  indicates  a  central-difference  operator.  The  U  and 
V  terms  are  the  velocities  described  previously  and  are 
found  as  shown: 

Oj+i/tA  *  {A+t  +  B+n  +  D)j+ i/tj, 

V},k+i/i  *  (B4t  +  +  E)j,k+  i/j 

where  the  A,B,C,D,  and  E  metrics  are  found  at  the 
node  points,  using  standard  second-order-accurate  finite- 
difference  operators.  All  metrics  quantities  are  found  and 
stored  at  node  points  (jjt),  and  simple  averages  are  used 
for  values  away  from  grid  nodes. 

The  density  values  used  in  the  continuity  equation, 
#>+>/>,*  PjM- i/s.  *»  defined 

+»/M  *  K1  ”  •'Mf+t/s,*  +  tyft/M  Py+ i+i/*,a 
PiA+i/i  - 1(1  -  v)p\i,k+i/i  +  •'y.s+i/a  Py.s+m+i/a 
where 

l  =*  ±1  when  $0 

m  *  ±1  when  Vy+i/M  >° 

Density  k  found  after  a  sweep  of  the  flow  field  has 
updated  all  the  A  values  from  the  nondimensional  equation 

p  -  [i  -  o.2m;(u*  +  k3-  r1)] *•*  (io) 

where  Af,  k  the  Mach  number  of  the  cylinder  surface. 
Borrowing  from  finite-volume  techniques,  the  densities  are 
found  at  cell-centered  locations  (j+1/2,  k+1/2),  which  are 
then  averaged  appropriately  to  yield  the  needed  values  at 
(i+l/2Jk)  and  at  (jJk+1/2).  Thk  technique  has  been  shown 
by  South  in  Safes  at  el.3*  to  yield  a  more  compact  calcu¬ 
lation  cell  for  density.  The  cell-centered  values  require  A 
derivatives  found  separately  from  the  derivatives  used  in 
the  continuity  aquation. 

The  switching  function  v  determines  when  to  start  us¬ 
ing  an  upstream  density  value  and  k  hence  controlled  by 
Mach  number.  Control  for  1/  in  the  (-direction  k  given  by 

‘'y+i/s  Jkm  "»«*  [0,1-1  /«£*]  for  0y+ i/m>°  (11) 
and 

•'y+i/M  “  me*[0,l-l/Af/+lt*]  for  l/y+i/*,*  <  0  (12) 

Control  in  the  vertical  direction,  or  q-direction,  k  given 
similarly. 

A  SLOR  technique  solves  the  linearised  continuity 
equation.  Linearisation  comes  from  simply  lagging  the  den¬ 
sity  one  iteration  step  so  that  for  each  iteration,  density  k 
considered  constant  and  k  updated  after  each  horisontal 
sweep  of  the  flow  field. 

Steger  and  Caradonna1*  show  an  alternative  to  the 
standard  switched  schemas  which  proves  stable  in  almost 
all  regimes: 


+  (12) 

In  thk  and  subsequent  sections  the  Vm  and  A.  symbok 
represent  backward  and  forward  difference  operators,  re¬ 
spectively,  with  the  subscript  denoting  direction.  Combi¬ 
nations  of  these  two  symbok  shew  higher  derivatives  such  es 
the  second-order  differences:  V,V,  for  a  backward  differ¬ 
ence  or  AaVa  for  a  central  difference.  Fortunately  thk  mix 
of  differencing  operators  closely  corresponds  to  the  artificial 
density  method  and  permits  an  easy  analysk  of  stability. 


I 


u  will  b«  shown  briefly.  Details  of  this  correspondence  in 
explained  in  ths  pnpor  by  St«f«r  sad  Caradonna.1* 


Stabilit 


sis  of  ths  Solution  Schtms 


A  stability  analysis  of  any  flnito-dilfsrsncs  msthod 
comas  Grom  ths  Von  Neumann  linear  stability  analysis 
ILSA).  Solving  for  ths  amplification  factor  rasulta  in  ths 
following  expression,  whose  magnitude  most  be  lees  than 
one  for  stability: 


{eo*(a)+t[ejna] }  / 

{2  -  coea  +  A/1  (I  +  co*2a  -  2 eo*a)+ 

4  .  (14) 

n  2 

•(etna  +  Jif*(2 ana  -  «tn2a)|) 

where  1  is  the  computational  cell  aspect  ratio  (Ay/Ax)1, 
and  a  and  P  are  the  wave  numbers  for  the  *-  and  in¬ 
directions,  respectively. 

The  usual  error-amplification  portraits  are  shown  in 
Fig.  4  for  the  undamped  and  damped  cases.  As  Mach  num¬ 
ber  increases,  an  aaeoriated  increase  in  amplitude  for  lower 
wave  numbers  ( low-frequency  errors)  is  seen.  The  stability 
limit  at  Mach  1.4  is  seen  by  the  flat  line,  showing  unity 
gain  for  errors  with  wave  numbers  near  sera.  Any  higher 
speeds  will  quickly  go  unstable  as  shown  in  the  curve  mm 
the  Mach  2.0  case. 

A  damping  term  added  to  the  iteration  or  correction 
matrix  allows  stable  solutions  with  higher  Mach  numbers. 
This  damping  term  also  increases  the  diagonal  dominance 
of  the  correction  matrix;  it  simply  adds  a  spatial  operator 
to  the  correction  term  in  the  general  relaxation  equation 
below: 

(S  +  DV.)C*  +  -  0  (IS) 

where  C*  is  the  correction  Ua+>  -  4");  ■  the  resid¬ 

ual,  which  indicates  how  well  the  finite-difference  equations 
have  converged;  and  w  is  the  overrelaxation  control,  which 
is  usually  set  within  the  range  1.9  •  2.0.  The  D  coefficient 
controls  the  amount  of  damping  added  to  the  solution.  This 
approach  will  not  change  the  solution  since  the  correction 
term  goes  to  sera  as  the  solution  converges.  Only  the  path 
of  the  convergence  changes,  so  that  stability  is  maintained. 
Hence,  a  new  time-path  is  taken  if  the  iteration  levels  are 
considered  time- levels.  This  variation  may  have  serious  im¬ 
plications  for  time-accurate  codes. 

Parametric  studies  of  this  damping  D  show  that  it 
must  increase  as  flf#/1  to  maintain  stable  iterations.  Tor 
the  cases  shown,  at  Mach  numbers  of  2.0  and  4 JO,  damping 
coefficients  of  -1  and  -4,  respectively,  maintain  stability. 
The  damping  effect  pulls  the  amplification  factor  curve  to 
the  unity  line  in  arose  where  it  was  previously  very  high. 

GridCfauntioa 


A  sheared  “H"  grid,  similar  to  grids  usually  used  for 
nonlifting  airfoil  studios ,  is  used  here  for  all  rectilinear  mo¬ 
tion  cases.  The  cylindrical  motion  cases  require  a  more 
sophisticated  grid  concentration  in  regions  where  the  shock 
pulse  is  expected  to  more  into  the  flow  field.  The  follow¬ 
ing  Poisson- like  equation  specifies  grid  concentration  in  one 
direction  only: 

-  C  F{t)  r 

The  C  factor  determines  how  much  the  grid  clusters  to  the 
control  function,  F(4).  The  control  function  varies  in  the 
4-direction,  so  that  the  grid  dusters  to  any  curve.  This 
curve  usually  comas  from  the  linear  characteristic  line  so 
that  grid  cells  cluster  to  the  region  where  Information  from 
the  airfoil  travels.  Ths  radial  lines  are  specified  separately. 
This  son-orthogonal  clustering  to  a  general  line  imparts 
a  skewing  of  the  grid  cells  that  creates  convergence  and 


similar  to  grids  usually  ussd  for 
s  used  here  for  all  rectilinear  mo- 


stability  problems  for  the  solution  algorithm.  Hence,  only 
a  weak  clustering  is  used. 

Boundary  Conditions 

Flow  tangescy  at  the  body  surface  is  enforced  by 

delta*  ■  “  C**)  I*- » 

The  far  boundary  conditions  specify  nonreflection  condi¬ 
tions  is  a  similar  to  those  used  fat  airfoil  codes.17 

Ths  general  potential  equation  is  solved  at  ths  outer  bound¬ 
ary  of  the  usual  arbitrary  and,  hence,  wrong  bound¬ 

ary  condition  of  setting  4  to  saro.  This  nonreflection  ap¬ 
proach  solves  an  approximation  shown  below  to  the  poten¬ 
tial  equation: 

0pd#  +  ?*r»O  (1®) 


where  0  and  V  are  const* 
ation;  they  are  defined  as 


found  from  the  previous  iter- 


0  m  Y+  — Vt4  +  U' 

V  "  (IT) 

Vrd  +  V' 

Ths  ( )'  quantities  change,  depending  on  the  boundary  lo¬ 
cation  as  defiwsd  by  tbs  following: 

U'  »  0,  V  »  -i-;  on  outer  boundary 


}  v,  } 

>/2V.’  y/iht. 


at  upper  roar  corner  point 


IT  «  V' m  0;  along  aft  boundary 
«• 

The  upstream  4  values  ate  held  to  saro,  and,  of  course, 
the  to  computational  variables  transforms  all  of  ths 

above  boundary-condition  equations  to  (,  q  space. 

The  nonrefiection  boundary  condition  above  [Eq.  (18)1 
produces  much  better  solutions  whan  proper  (upetream) 


produces  much  better  solutions  whan  proper  (upetream) 
difference  operators  are  used  as  naadad.  This  upstream 
differencing  is  sufficiently  implemented  by  using  backward- 
difference  operators  on  both  the  vertical  and  horiaontal 
derivatives  and  using  a  marching  procedure  toward  the  up¬ 
per  right-hand  corner  point. 


The  computational  code  was  first  validated  in  the 
cartesian  mode  and  compared  with  a  similar  code's  results. 
The  validated  code  was  than  used  in  the  cylindrical  mode 
to  predict  the  experimental  results.  After  finding  differ¬ 
ences  in  the  experimental  results  and  the  code  predictions, 
corrections  to  the  code  were  used  to  modify  the  predicted 
results. 

Computational  Rasulta 


This  cods  shews  proper  solutions  when  compared  with 
other  established  transonic  codes.  Rectilinear  motion  re¬ 
sults  compare  well  with  similar  results  published  by  Holst 
and  Ballhaus11  for  Mach  numbers  of  0.70  and  0.84  using 
a  circular-arc,  10%- thick  airfoil.  Stability  limitations  are 
a  problem  for  potential  codas  like  the  code  of  this  study 
when  the  flow  field  contains  a  region  of  high  speed  flow 
(e.g.,  Mach  1.4  or  greater).  The  previous  section  on  sta¬ 
bility  shows  how  a  linear  model  of  the  potential  equation 
stabilises  at  Mach  numbers  of  2.0  and  greater  for  properly 
controlled  damping.  A  high- Mach- number  case,  which  is 
stabilised  by  this  added  damping,  appears  next.  Figure  S 
shows  the  surface  pressures  from  a  Mach  2.0  rectilinear 


motion  test.  The  solid  Um  in  Pig.  S  shows  Um  txptcUd 
linaar  results  for  this  supersonic  caaa  aa  (bond  from  shock- 
nptmkiB  theory.  The  analytical  and  computed  aointiona 
match  vary  trail,  thereby  demonstrating  that  difficult  and 
usually  unstable  cases  are  poeaibie. 

The  results  in  Tig.  8  come  from  a  high- Mach- number 
case  to  shots  how  added  damping  controls  an  inherent  in* 
stability.  The  same  instability  problems  arise  in  the  cylin¬ 
drical  casea  when  the  velocity  given  by  Or  becomes  greatly 
supersonic.  Damping  parameter  levels  far  thia  Mach  2.0 
case  are  about  what  are  expected  from  the  linear  model 
case,  bat  only  near  convergence.  Daring  the  convergence 
process,  the  solution  error  starts  to  grow  quickly,  requiring 
an  increase  in  the  damping  coefficient.  Automatic  damping 
adjustment  avoids  guesswork  in  picking  a  damping  value. 
This  adjustment  algorithm  is  shown  if  D  is  taken  as  the 
damping  value  and  superscripts  indicate  iteration  count: 

D“  -1.06  D—‘  if  Rm  >  1.00/T-1  ,  . 

(18) 

D*  —0.90  D“"‘  if  R*  <  OMR*'1 

The  damping  magnitude  is  usually  constrained  to  a  value 
of  at  least  0.3  but  lees  than  25. 

Complete  delocalisation  occurs  for  a  Mach  number  of 
0.79  as  shown  in  Fig.  6b  by  the  radiating  wave.  Parameters 
for  this  case  consist  of  an  aspect  ratio  of  7,  an  airfoil  thick¬ 
ness  (r)  of  129f,  and  surface  Mach  numbers  of  about  0.75. 
These  physical  parameters  match  those  from  the  experi¬ 
mental  apparatus.  These  results  are  not  unlihs  those  seen 
in  the  coane-grid  solutions.  The  fine-grid  surface  pressures 
seen  in  Fig.  6a  show  a  finer  shock  width  and  a  slightly 
higher  negative  peak  than  the  coarse  grid  results.  The  de¬ 
localised  wave  seen  in  Fig.  fib  is  aiao  resolved  much  better 
in  this  fine-grid  case. 

The  shock  on  the  surface  still  damps  out  to  farm  a 
fairly  smooth  pulse  away  from  the  airfoil.  This  damping  is 
not  entirely  due  to  a  loea  of  resolution  from  grid  coarsen  aw. 
Grids  of  similar  node  points  used  on  narrower  domains  yield 
a  finer  mesh  but  do  not  appreciably  improve  the  final  so¬ 
lution.  A  grid  spacing  of  20  points  par  chord  length  is 
believed  to  be  sufficient.  Covering  tha  flow  field  out  to  the 
16-chords  location  with  20  points  par  chord  length  in  two 
directions  uses  almost  all  of  the  conventional  memory  of 
2  x  10*  words  on  the  Cray  XMP-4/8  at  the  Ames  Research 
Center.  A  finer  grid  would  require  lengthy  run  times  and 
special  coding  procedures  outside  tha  scope  of  this  work. 

Reralling  that  the  wave  equation  in  two-dimensions 
can  not  support  a  sharp  discontinuity,  such  as  a  shock 
explains  why  the  shock  should  damp  out  aa  it  moves  off 
the  surface.  The  three-dimensional  wave  equation  does 
support  and  propagate  impulses  though.**  Therefore,  tha 
sharp  pulses  seen  in  hovering  rotor  testa  are  not  expected 
in  the  code  results  or  in  the  two-dimensional  experimental 
results. 

The  wavelet  that  exits  at  the  rear  boundary  shows 
how  well  the  nonreflection  conditions  imposed  at  the  outer 
boundaries  work  to  accurately  handle  the  flow  near  a 
boundary.  Canes  that  do  not  have  this  condition,  but  In¬ 
stead  use  the  usual  procedure  of  simply  specifying  the  4 
value,  produce  many  reflected  waves  at  the  boundaries  far 
even  simple  cases.  A  clearly  evident  wave  reflection  ap¬ 
pears  when  such  arbitrary  boundary  conditions  are  applied 
and  when  the  outer  boundary  grid  is  sufficiently  fine  to  re¬ 
solve  the  reflected  wave.  Another  Indication  of  incorrect 
boundary  conditions  is  seen  by  a  build-up  of  4  values  in 
the  upper  right-hand  corner  of  the  solution  field.  This 
effect  occurs  in  three  dimensions  also.  When  the  nonre¬ 
flection  condition  is  disabled  in  another  fall-potential  but 
thiee-dimeasianai  code,  a  similar  accumulation  occurs  at 
the  outer  boundary.**  The  abeorbtion  boundary  condition 
produces  remarkably  improved  solutions  near  the  bound¬ 
aries  by  eliminating  both  of  these  problems.  The  reflection 


error  is  not  a  problem  far  most  codes,  because  the  grid  near 
the  outer  boundaries  is  usually  sufficiently  coarse  to  damp 
out  erroneous  boundary  reflections  that  ndght  cause  signifi¬ 
cant  error  in  the  surface  solution,  and  the  outer  boundaries 
are  usually  placed  very  far  from  the  body. 

The  results  found  from  the  fine-grid  solution  presented 
in  Fig.  6  are  reformulated  to  produce  the  results  shown  in 
Fig.  7.  This  presentation  corresponds  to  the  results  ex¬ 
pected  from  aricrophone  date  taken  in  terms  of  relative 
pressure  during  the  experiment.  The  microphones  were  lo¬ 
cated  0.5, 1.0, 34),  and  6.0  chorda  away  from  the  surface. 

Experimental  Results 

The  two-dimensional  experiment  described  hare  repro¬ 
duces  Boost  of  ths  conditions  awnmsrl  fay  the  computer 
code.  A  boundary-layer  correction  to  the  code  allows  for 
tha  viscous  effects  seen  in  tha  experiment.  A  preliminary 
computer  pjamltetoi  of  this  tote  showod  that  doiocoUsod 
flow  would  occur  at  Mach  0.81,  dictating  an  operating  speed 
of  4973  rpm.  Although  the  final  simulations  show  delocal¬ 
isation  at  Mach  0.79,  the  higher  speed  was  used  in  the 
design  of  the  cylinder.  Figure  3  shows  a  sketch  of  the  back 
setup  of  a  cylinder  spinning  between  two  walls  to  maintain 
a  two-dimensional  flow.  Pressure  sensors  wars  mounted  in 
one  support  wall  along  a  radial  Uns  from  the  cylinder  origin, 
and  the  sensor  faces  ware  kept  flush  with  the  wall  surface. 
The  preliminary  design  of  the  cylinder  was  done  at  Asm 
Research  Canter.*1 

Figure  8  shows  1)  ths  field  pressures  sway  from  the 
rotating  cylinder  far  part  of  a  cylinder  revolution  and  2) 
how  the  pulse  magnitude  decreases  with  radial  distance. 
All  of  tha  plots  from  experimental  data  come  from  a  teat  in 
which  the  surface  Mach  number  was  0.824.  This  velocity  is 
wall  above  ths  predicted  delocalisation  speed  yet  not  a  hint 
of  delocalisation  is  aaan.  A  pressure  pulse  similar  to  that 
expected  from  the  computer  code  at  Mach  0.79  is  asm,  but 
the  magnitude  is  low  by  a  factor  of  5.  This  reduced  pulse 
magnitude  is  due  to  the  thick  boundary  layer  observed  in 
the  test  despite  all  efforts  to  reduce  viscous  affects  through 
the  boundary- layer  removal  system.  Since  ths  boundary- 
layer  reduces  the  velocity  of  the  body  relative  to  the  moving 
air  just  above  the  body,  leas  of  a  net  disturbance  is  created. 

Figure  9  directly  comperes  experimental  pressure  re¬ 
sults  at  the  first  transducer  location  of  04  chord  far  a  fall 
revotution  with  ths  single  pulss  from  ths  computer  predic¬ 
tion.  Tha  observed  pulse  shape  is  about  the  same  as  ex¬ 
pected  from  the  pulse  prediction.  No  strong  shock  is  seen 
in  either  result,  although  the  pulse  is  not  symmetrical  since 
it  shows  a  sharper  recovery  to  free  atrcom  pressure.  The 
smaller  spike  and  large  area  of  slightly  positive  pressure 
just  forward  of  the  main  negative  pulse  in  the  test  data 
comas  from  variations  from  the  circular  that  were  present 
in  the  cylinder.  Clearly,  the  code  must  include  both  the 
exact  cylinder  shape  and  the  boundary-layer  influences. 


Boundary  -  Layer  Corrections 


Tha  approximate  boundary- layer  correction  added  to 
the  inviscid  potential  code  comm  from  aa  analytical  solu¬ 
tion  to  the  Navier-Stokm  equations.  This  correction  in¬ 
volves  modifying  the  fir  component  of  the  rotating  coor¬ 
dinates  transformation  so  that  the  no-elip  condition  is  met 
at  the  surface  and  an  appropriate  velocity  profile  blends  to 
the  inviscid  free  stream  Ths  assumptions  in  this  deriva¬ 
tion  are  incompressibility,  laminar  and  steady-state  flow, 
and  uniform  suction  instead  of  the  discrete  suction  ports. 
An  eddy-viscosity  model  takas  care  of  the  turbulence  mod¬ 
eling. 

Allowing  suction  and  looking  for  a  steady-state  solu¬ 
tion  with  v  ■  RV/r  from  the  continuity  equation,  where  R 
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ia  the  cylinder  radius  and  V  ia  the  uniform  auction  velocity 
at  the  surface,  producea 


where  a  ia  the  velocity  component  in  the  ^-direction, 
and  v  ia  the  kinematic  viacoaity  coefficient.  The  auction 
Reynoida  number  ia  defined  aa  R,  *  VRIv,  which  ia  nega- 
tive  becauae  the  auction  velocity  ia  negative,  and  the  cylin¬ 
der  angular  apeed  ia  given  by  fl.  Various  teats  show  that 
the  pressures  can  be  well  predicted  by  using  R,  a  —4. 
This  value  yields  the  same  pressure  (-1.1  lb/in.a)  that  was 
measured  in  the  experiment  at  the  0.5-chord  location  and 
producea  a  reasonable  looking  flow-velocity  profile.  This 
modified  velocity  introduces  a  minor  variation  from  the  as¬ 
sumed  irrotational  flow  field,  which  predicts  density  with  a 
negligible  5/5  error.  However,  this  correction  only  matches 
the  magnitude  of  the  main  pulse  and  fails  to  match  to  over¬ 
all  pulse  shape  seen  in  Fig.  9. 

A  detailed  model  of  the  geometric  boundary  condi¬ 
tion  on  the  cylinder  surface  reflects  the  exact  variations 
from  circular.  Pressure  results  using  the  exact  geometric 
conditions  for  the  cylinder  are  shown  in  Fig.  10  for  a  do¬ 
main  of  a  little  over  one-half  of  a  revolution.  The  pressure 
results  now  show  a  great  deal  of  similarity  to  measured  ex¬ 
perimental  pressures  (Fig.  11),  except  that  the  computed 
magnitudes  are  still  high  by  a  factor  of  5.  The  smaller  pulse 
after  the  large  pulse  from  the  airfoil  bump  is  now  well  mod¬ 
eled,  and  the  slowly  rising  positive-pressure  hump  observed 
in  the  test  is  now  seen.  Mice,  an  accurate  model  has  been 
developed  to  predict  the  form  of  the  pressures  seen  dur¬ 
ing  the  experimental  tasting.  Now  only  a  correction  to  the 
magnitudes  is  required. 

This  exact  cylinder  shape  case  uses  the  same  R»  «  -4 
exponent  that  waa  previously  found  to  produce  such  good 
matches  to  the  observed  pressures.  Complete  details  at  how 
the  potential  code  includes  these  boundary-hirer  influences, 
as  well  as  the  exact  boundary  conditions  o f  the  cylinder, 
are  given  in  Ref.  30.  The  results  from  such  a  fully  cor¬ 
rected  case  are  shown  in  Fig.  12  where  they  are  compared 
with  the  observed  experimental  results  for  the  0.5-chord  lo¬ 
cation.  The  boundary  layer  has  smoothed  over  the  sharp 
peaks  seen  in  the  in  viscid  prediction.  However,  the  pulse 
magnitudes  are  much  closer  than  those  seen  ia  Fig.  11.  The 
addition  of  exact  boundary  conditions  increases  the  main- 
peak  magnitude,  but  there  has  been  no  attempt  to  correct 
this  peak.  Since  surface  details  are  most  important  in  thin 
case,  a  grid  is  used  that  is  clustered  to  the  body  with  aa 
almost  uniform  grid  in  the  circumferential  direction.  The 
computational  results  at  points  farther  away  from  the  sur¬ 
face  show  a  definite  loss  of  resolution,  although  the  results 
farther  out  show  the  same  trends  of  smeared  pulses  seen  ia 
the  0.5-chord  plot  of  Fig.  12. 


This  section  derives  the  observed  two-dimensional  re¬ 
sults  so  that  they  compare  with  expected  results  from  a 
real  helicopter  rotor.  Figure  13  shows  results  from  this 
two-dimensional  code  in  terms  of  helicopter  rotor  nomen¬ 
clature.  The  two-dimensional  “aspect  ratio”  used  here  sink' 
ply  means  the  radius  of  the  rotating  cylinder  divided  by  the 
chord  length  of  the  body  attached  to  the  cylinder’s  edge. 
This  aspect  ratio  is  similar  to  the  usual  helicopter  definition 
of  rotor  radius  divided  by  blade  chord. 

The  general  effect  of  various  aspect  ratios  and  thick¬ 
ness  ratios  is  not  surprising:  to  obtain  a  higher  deloralisa- 
tiou  Mach  number,  simply  disturb  tbs  air  Ism.  This  basic 
principle  can  be  achieved  by  using  a  thinner  body  or  a 
large  two-dimensional  aspect  ratio,  which  cones  ponds  to  a 
long,  slender  rotor  blade  on  a  real  helicopter.  Using  thinner 
bodies  is  rather  obvious  since  a  thhmor  body  is  biown  to 


produce  a  smaller  disturbance.  A  large  aspect  ratio  means 
that  the  relatively  short  body  must  move  faster  to  disturb 
the  air  enough  to  send  a  disturbance  off  the  tip  to  the  sonic 
circle.  Additionally,  larger  aspect  ratios  also  mean  there  is  a 
greater  distance  from  the  blade  tip  to  the  linear  souk  circle 
when  measured  in  chords.  The  linear  sonic  circle  is  defined 
as  a  line  where  the  fir  term  of  the  rotational  transforma¬ 
tion  equals  the  free  stresm  speed  of  sound.  The  distance 
from  the  cylinder  edge  to  this  sonic  circle  is  given  in  chords 
by  AL(1/Mt  -  1)<  which  increases  with  aspect  ratio  for  a 
given  cylinder  surface  Mach  number  A#«.  A  larger  aspect 
ratio  requires  a  larger  disturbance  to  bridge  the  increased 
gap. 

Figure  13  shows  a  consistent  dependence  on  aspect  rar> 
tio  for  most  airfoil  *  Sicknesses.  The  body  thirknooess  typ¬ 
ical  of  most  helicopter  blades  produce  curvee  that  are  al¬ 
most  parallel.  The  very  thin  bodice-395  and  lam  show  a 
lack  of  dependence  on  aspect  ratio.  Since  the  edge  pertur¬ 
bation  is  so  small,  connection  to  the  sonk  circle  requires  a 
very  high  Mach  number.  This  high  Mach  number  throws 
disturbances  off  the  body  that  quickly  develop  into  Mach 
waves,  a  result  of  nonlinear  acoustk  velocity  effects.  These 
Mach  waves  connect  to  the  sonk  circle  to  delocalise  the 
flow.  Since  delocalisation  comm  from  such  a  Mach  waive, 
the  distance  in  chords  is  irrelevant.  Small  disturbances  of 
various  lengths  produce  almost  identical  waves.  With  this 
lorn  of  chord-length  dependence,  the  distances  in  rotor  radii 
become  important.  Them  radii  distances  ate  given  m  a 
function  of  cylinder  Mach  number  jJ-  — 1.  The  delocalisa- 
tion  Mach  number  now  depends  only  on  this  distance  and 
hence  on  cylinder  Mach  number,  which  explains  the  almost 
flat  curve  aa  aspect  ratio  changes  for  very  thin  bodies. 

Since  the  above  approach  is  strictly  two-dimensional, 
tip  effects  are  not  included.  A  comparison  with  three- 
dimensional  results  reveals  when  tip  effects  are  important 
and  what  effects  the  tip  relief  will  have.  There  are  very 
few  test  date  available  on  delocalisation  Mach  number 
for  a  broad  range  of  rotor  parameters.  Hence,  a  small- 
disturbance  rotor  code  supplies  a  three-dimensional  date 
bom  ia  the  place  of  experimental  evidence.*1  Figure  14 
shows  the  results  from  this  three  dimensions!  code  tor  a  set 
of  parameters  similar  to  thorn  that  are  reflected  in  Fig.  13. 
Both  of  the  studies  mentioned  in  this  section  use  a  symmet¬ 
ric  parabolk-arc  airfoil  section.  The  full  rotor  prediction 
shows  very  similar  trends  to  thorn  from  the  two-dimensional 
approach  seen  in  Fig.  13.  The  results  for  the  thicker  blade 
show  a  smooth  dependence  on  blade  thickness  and  aspect 
ratio.  The  average  Mach  number  is  higher  for  the  three- 
dimensional  case,  out  the  thicker  blades  show  the  same  con¬ 
vex  curve.  As  the  blade  gets  thinner,  the  same  flattening 
of  the  curve  is  seen.  The  695- thick  blade  result  is  slightly 
twmw  to  tint  Um  mtiwhar  actnilljr 

increases  for  smaller  aspect  ratios.  The  two-dimensional 
remits  do  not  show  this  effect  until  blades  that  are  thinner 
than  295  are  used. 

The  tip  relief  in  the  three-dimensional  cases  explains 
both  the  higher  overall  delocalisation  Mach  number  and 
the  delay  in  curve  shape  reversal.  The  extra  dimension 
present  in  the  full  rotor  cases  allows  extra  room  for  dis¬ 
turbances  to  dissipate.  Hence,  a  higher  Mach  number  is 
required  in  the  three-dimensional  casss  to  achieve  a  distur¬ 
bance  strong  enough  to  delocalise  the  flow.  Since  higher  de- 
localisation  Mach  numbers  appear  overall,  the  Mach  wave 
delocalisation  discussed  above  becomes  important  for  the 
thkhar  blade  results,  such  aa  the  695  one  which  is  concave. 
Again  the  change  in  curve  shape  occurs  at  Mach  numbers 
of  0.9. 

Conclusions 

Techniques  were  developed  for  application  to  potential- 
flow  computational  fluid  dynamics  methods  for  predkting 


acoustic  wars*.  Using  the  potential-flaw  equation  to  define 
the  outer  boundary  values  for  the  potential,  d,  alleviates 
one  problem  of  reflected  waves.  A  stability  limit  problem 
occurs  at  Mach  numbers  above  1.4,  but  adding  a  damping 
factor  to  the  relaxation  equation  allows  higher  speeds.  This 
damping  factor  comes  from  a  linear  equation  analysis,  and  a 
demonstration  of  how  it  controls  divergence  was  presented. 

The  computational  results  show  a  definite  change  in 
the  flow  field  when  delocalisation  occurs.  The  pressures 
found  away  from  the  spinning  body  do  not  show  the  shock¬ 
like  waveform  seen  in  hovering  model-rotor  tests,  since  a 
two-dimensional  wave  solution  tends  to  damp  discontinu¬ 
ities  as  opposed  to  the  three-dimensional  solution  which 
preserves  spikes  in  pressure. 

A  two-dimensional  experiment  was  conducted,  and 
limited  experimental  results  were  presented  to  support  the 
findings  of  the  computer  prediction.  A  set  of  computer-code 
corrections  for  the  exact  cylinder  shape  and  the  boundary- 
layer  entrainment  produces  results  similar  to  those  seen  in 
the  experiment. 

Parametric  variations  with  the  two-dimensional  com¬ 
puter  code  show  a  smooth  dependence  on  airfoil  thirlmass 
and  aspect  ratio  for  configurations  similar  to  actual  heli¬ 
copter  rotors.  An  unusual  lack  of  aspect-ratio  influence 
occurs  for  very  thin  airfoil  sections  (e.g.,  395  and  less)  for 
which  thinner  disturbances  create  delocalisation  just  by 
starting  a  Mach  wave  that  is  independent  of  chord  length 
and,  hence,  aspect  ratio  fradius/cbord). 

Higher  delocalisation  speeds  overall  are  seen  in  the 
three-dimensional  rotor  code  owing  to  tip-relief  effects, 
which  reduce  any  tip  disturbance  and  thereby  require  a 
higher  speed  to  delocalise  the  flow.  The  reversal  of  de¬ 
pendence  on  aspect  ratio  occurs  for  thicker  blades  in  the 
rotor  case  but  for  about  the  same  Mach  number  as  the 
two-dimensional  model  case.  This  indicates  that  at  such  a 
high  speed,  Mach  waves  are  sufficient  to  delocalise  the  flow. 

The  model  results  and  the  rotor  results  show  a  clear 
correlation.  The  region  where  tip-relief  effecta  are  impor¬ 
tant  is  shown  and  the  two-dimensional  model  is  shown  to 
predict  the  expected  delocalisation  Mach  number  for  a  real 
helicopter  rotor  baaed  on  the  three-dimensional  rotor  code 
correlations. 
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